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Located in a dense complex environment, W51C provides an excellent scenario to probe accelerated protons in 
SNRs and their interaction with surrounding target material. Here we report the observation of extended Very 
High Energy (VHE) gamma-ray emission from the W51C supernova remnant (SNR) with MAGIC. Detections 
of extended gamma-ray emission in the same region have already been reported by the Fermi and H.E.S.S. 
collaborations. Fermi/L AT measured the source spectrum in the energy range between 0.2 and 50 GeV, which 
was found to be well fit by a hadronic neutral-pion decay model. The VHE observations presented here, 
obtained with the improved MAGIC stereo system, allow us to pinpoint the VHE gamma-ray emission in the 
dense shocked molecular cloud surrounding the remnant shell. The MAGIC data also allow us to measure, 
for the first time, the VHE emission spectrum of W51C from the highest Fermi/LAT energies up to TeV. 
The spatial distribution and spectral properties of the VHE emission suggest a hadronic origin of the observed 
gamma rays. Therefore W51C is a prime candidate for a cosmic ray accelerator. 



1. Introduction 

Supernova Remnants (SNRs) are believed to be the 
sources of galactic cosmic rays for many years. The 
underlying process from a theoretical point of view is 
the diffusive shock acceleration in which charged par- 
ticles are accelerated in the expanding shocks of the 
SNRs, receiving part of the kinetic energy of the shock 
through the first order Fermi acceleration [2]. Even 
though the last years provided increasing support by 
observational evidences, there is still no absolute proof 
that SNRs can accelerate hadronic particles. The ob- 
servations of high energy (HE, > 100 MeV) and very 
high energy (VHE, > 100 GeV) gamma rays from 
many SNRs clearly showed that they are sites of rela- 
tivistic particles of energies compatible to the emission 
seen in gamma rays. To unambiguously establish the 
acceleration of hadronic particles, the hadronic ori- 
gin of the observed gamma rays needs to be proven. 
In almost all cases, electromagnetic scenarios (gamma 
rays produced by up-scattering of photons through 
Inverse Compton scattering by accelerated electrons) 
and hadronic scenarios (gamma rays produced after 
the decay of 7r°) cannot be distinguished by looking 
only at the VHE gamma-ray regime. Therefore obser- 
vations at other wavelengths are needed to better un- 
derstand the production mechanism of gamma rays. 
This includes on the one hand measurements of the 
parameters describing the interstellar medium (ISM) 
around the object like magnetic field and density. On 
the other hand observations of HE gamma rays in the 
range from MeVs to GeVs are very important because 
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they can reveal features in the spectrum that can dis- 
tinguish both possibilities. Since the luminosity of 
gamma ray emission of a hadronic origin is propor- 
tional to the product of CR density and the density of 
the surrounding ISM, the existence of dense molecular 
clouds close to the SNR acts as a natural amplifier for 
the hadronic component. A tool to detect evidence 
of cosmic rays acceleration is therefore provided by 
gamma rays emitted by molecular clouds adjacent to 
or shocked by SNRs Q. W28 and IC443 [3 are 
two examples where this effect might be responsible 
for the production of VHE gamma rays. In the W28 
case, gamma rays might be produced by cosmic rays 
that escaped the accelerating SNR while in IC 443 a 
Crushed Cloud scenario might be at work where the 
molecular cloud is shocked by the blast wave region 
of the SNR [6|. A major differ nee in the two sce- 
narios is the fact that the illumination of molecular 
clouds by runaway cosmic rays is affected by the en- 
ergy dependent escape and diffusion of cosmic rays 
and might therefore show spectral features due to this. 
In the case of a Crushed Cloud the diffusion effects are 
thought to be much smaller and the emission could 
be produced by the cosmic ray source spectrum itself. 
On the other hand re-acceleration of cosmic rays in 
the shocked cloud may play a role in this case. 

The SNR W51C is a promising candidate to test 
and study cosmic ray acceleration in SNRs. W51C 
is a mixed morphology SNR with an elliptical shape 
in radio and a size of 0.8° x 0.6°. It is located in 
the tangential point of the Sagittarius arm at a dis- 
tance of ~ 6 kpc [i[ and the estimated age is around 
30 kyrs. The W51 complex consists of three main 
components: Two star- forming regions, W51A and 
W51B, and the SNR W51C. While W51A is separated 
from the other two, W51B overlaps with the North- 
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Western rim of W51C. Shocked atomic and molecular 
gases have been observed in radio data [Toj| [llj ]. pro- 
viding direct evidence on the interaction of the W51C 
shock with a large molecular cloud. This strongly indi- 
cates that W51C represents a Crushed Cloud scenario 
rather than a cloud illumination scenario. In order 
to distinguish the two cases, a precise knowledge of 
the spatial origin of the VHE emission is essential. 
X-ray emission has been detected from the star form- 
ing regions W51A and W51B. W51C is also visible 
in X-rays showing both a shell type and center-filled 
morphology [7[. A non-thermal X-ray emission has 
only been detected from the relatively bright and com- 
pact source CXO J192318. 5+140505, which is thought 
to be a pulsar wind nebula (PWN) associated to the 
SNR H. 

Fermi/LAT detected gamma-ray emission from 
200MeV to 50GeV extended throughout W51B and 
W51C [l2|. The relatively large point-spread- function 
PSF does not allow to tell which objects of the field 
of view is responsible for the observed emission. The 
Fermi/LAT emission region is extended in comparison 
to the PSF of the instrument (cr ex t = 0.22 ± 0.02 deg. 
The H.E.S.S. Collaboration detected VHE emission 
above 420 GeV from an extended region coincident 
with the Fermi/LAT emission region [13]. In their 
sky map, the H.E.S.S. emission is smoothed with a 
radius of 0.22° [l3| and, as well as in the case of 
Fermi/LAT, it overlaps with several HII regions, the 
molecular mass in W51B as well as the PWN can- 
didate CXOJ192318.5+140305. The flux measured 
by H.E.S.S. above 1 TeV is at the level of 3% of the 
Crab Nebula flux. Finally, also the MILAGRO Col- 
laboration reported a possible excess from the same 
source at energies above several TeV [14j|. The mod- 
eling of the spectral energy distribution (SED) mea- 
sured by Fermi/LAT suggests a hadronic mechanism 
as the main origin of the gamma rays. The interaction 
should take place in the region between the supernova 
remnant and the dense, complex region North- West of 
it. In this case, the engulfed cloud would be those de- 
scribed in [3 EH- 



2. Observations with the MAGIC 
telescopes 

MAGIC consists of two 17 m diameter imaging at- 
mospheric Cherenkov telescopes located at the Roque 
de los Muchachos in the Canary Island of La Palma 
(28°46 / N, 17°53'W) at the height of 2200 ma.s.l. As- 
tronomical observations of VHE gamma-ray sources 
are performed with the two telescopes simultaneously 
which provides a major improvement of performance 
with respect to the single telescope observations pre- 
viously done with MAGIC-I [16]. MAGIC observed 
W51C between May 17 and August 19 2010. The 
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Figure 1: Distribution of the squared angular distance 
between the reconstructed VHE gamma-like events (red 
points) above 150 GeV and the source position together 
with background (black points) reconstructed from a 
symmetric position in the camera(s). The vertical dotted 
line marks the region defined to compute the excess 
events. 



observations were carried out in the so-called wob- 
ble mode and covered a zenith angle range between 
14 and 35 degrees. The Fermi/LAT source W51C 
(RA=19.385 h, 5=14.19 °) was chosen as central po- 
sition for the observations. After applying quality cuts 
we collected a total of 31.1 h effective 1 dark time. All 
data were taken in stereoscopic mode, recording only 
events that triggered both telescopes. The trigger en- 
ergy threshold of the system is around 50 GeV [161 ] . 
This is the lowest energy threshold among IACTs in 
operation and provides the chance to have ground 
based observations with an energy range overlapping 
with the one of Fermi/LAT. The analysis of the data 
was performed using the MARS analysis framework 
which is the standard software used for MAGIC data 
analysis [l7j- The details of the analysis, as well as the 
general performance of MAGIC in stereoscopic mode, 
are reported in [l6| [l5| . 



3. Results from MAGIC observations 

MAGIC founds a spatially extended excess of VHE 
gamma-ray events from the direction of W51C. In to- 
tal 924 excess events were found above 150 GeV in 



-'^Due to dead time of the MAGIC-II telescope readout sys- 
tem the effective observation time is always lower than the real 
observation time. 
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Figure 2: VHE gamma-ray emission from W51C obtained 
with the MAGIC telescopes above 150 GeV. The map 
has been smoothed with a Gaussian kernel of cr=0.10°. 
The color scale shows the relative flux of gamma rays 
(excess events normalized to the number of background 
events) and the black contours are different levels of the 
test statistic variable (Li & Ma eq. 17 [l9[ applied on a 
smoothed and modeled background estimation. It 
roughly corresponds to a Gaussian significance. More 
details given in [18[). In the figure the region defined for 
integrating the signal (dotted line) and the PSF after 
smearing (bottom left inset) are also shown. 



the 31.1 hours of effective time (see figure [T]), result- 
ing an a detection significance of 8.29 a. The source 
as seen by MAGIC shows an extension (sigma of two- 
dimensional Gaussian fit) of 0.16±0.02°. This lies well 
above the PSF of the analysis (0.08° in the same en- 
ergy range assuming a spectral index of -2.6). The 
centroid of the emission detected by MAGIC is coin- 
cident with the centroid of the Fermi/LAT emission: 
RA=19.387±0.002 h, £=14.18±0.02°. A sky map of 
the MAGIC observation at energies above 150 GeV is 
shown in figure [2j Details about MAGIC sky maps 
and test statistic are given in [18|. 

MAGIC has measured the differential energy spec- 
trum of the VHE gamma-ray emission in the energy 
range 75 GeV to 3.3 TeV. The differential spectrum 
is well fitted by a power law as can be seen in figure [3] 
(x 2 /ndf = 4.5/5). The obtained spectral index of the 
VHE gamma-ray emission is -2.40±0.12 s tat- The in- 
tegrated flux above 1 TeV corresponds to 3.8% of the 
Crab Nebula in agreement with the integral flux re- 
ported by H.E.S.S. above 1 TeV. The obtained 
flux in units of TeV -1 cm -2 s _1 is given by: 
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Figure 3: Measured flux from W51C with the MAGIC 
telescopes. The red points show the differential fluxes. 
The red line represents the best fitted power-law. As a 
reference the measured by MAGIC-I differential energy 
spectrum from the Crab Nebula is also shown. 



4. Discussion 

The MAGIC data do not only fill the gap between 
the Fermi/LAT and the H.E.S.S. measurements but 
furthermore show that the spectral distribution can 
be described by a single power-law from energies of 
~10 GeV up to TeV energies. Figure |4] shows the SED 
measured by Fermi/LAT and MAGIC together with 
the H.E.S.S. measurement converted into a differential 
flux 2 . MAGIC data agrees well with the Fermi/LAT 
and H.E.S.S. measurements. In the same figure the 
predictions from the phenomenological model used by 
Fermi/LAT in [12| to explain the origin of the gamma- 
ray emission are also shown. Three different scenar- 
ios are considered: one where gamma-ray emission 
is dominated by 7r° decay and two more where the 
gamma-ray emission is mainly dominated by inverse 
Compton (IC) or Bremsstrahlung of electrons. 

The VHE MAGIC data points as well as H.E.S.S. 
point are above the predictions from the model in 
all scenarios at energies above 1 TeV. Although this 
seems to favor the electromagnetic scenarios, both lep- 
tonic models show difficulties in fitting the radio data 
and, in the case of the IC model, severe energy con- 
straints. In addition both leptonic models need an 
electron to proton ratio in order of one. Such a value 
is much higher than the one observed in direct cos- 
mic ray measurements (e/p ~ 10 -3 ) and difficult to 
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2 The integral flux above 1 TeV reported in [13| is converted 
into a differential flux using the MAGIC spectral index of -2.4. 
An error of ±0.4 is assumed in order to obtain the error bars 
shown. 
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5. Conclusions 
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Figure 4: SED of W51C measured by Fermi /LAT (black 
points), MAGIC (red points) and H.E.S.S. (green point). 
Also shown in the figure the 3 different scenarios for 
modeling the multi wavelength data that were used 
in [l2j]. Gamma-ray emission is explained by emission 
coming from a hadronic dominated scenario (continuous 
blue line), inverse Compton dominated scenario 
(dot-dashed line) or a Bremsstrahlung dominated 
scenario (dashed line). 



MAGIC observations do not only confirm the emis- 
sion of HE-VHE gamma rays from an extended source 
located in the SNR W51C, but extend the spectrum 
almost without gap above the Fermi /LAT energies. 
The measured SED is in agreement with a extrap- 
olation of the Fermi data and agrees with the flux 
measured by the H.E.S.S. collaboration. The emission 
measured by MAGIC is spatially coincident with that 
reported by Fermi/ 'LAT. Moreover, the higher angular 
resolution provided by MAGIC shows that the bulk of 
the VHE gamma-ray emission comes from the shocked 
molecular cloud located where the SNR shock engulfs 
a large molecular cloud, creating a shocked gas region 
distinguishable in the radio data. This fact, a constant 
spectral index of the emission from ~GeV to multi- 
TeV, and the problems to explain the emission in a 
leptonic scenario suggest that the gamma-ray emis- 
sion has most likely a hadronic origin as it is expected 
in the case of gamma rays produced in a crushed cloud 
scenario. 



explain due to the much higher energy losses for elec- 
trons. Concerning the hadronic emission the param- 
eters used by the Fermi /LAT- team leads to a proton 
spectral index between ~10 and ~500GeV close to 
-3. A slightly harder index would be able to explain 
both MAGIC and the Fermi /LAT points. The spec- 
tral index measured by Fermi /LAT above 10 GeV 
is —2.5 ±0.18 [l|. Taken this into account the non 
matching of this hadronic scenario does not rule out a 
hadronic origin of the emission. Moreover a hadronic 
origin predicts a constant slope of the gamma emission 
above several GeV up to the cut-off originating from 
reaching the maximum energy of protons responsible 
for the emission. 

The angular resolution of MAGIC at energies of 
100 GeV is comparable to that of Fermi/ 'LAT. For en- 
ergies above 700 GeV, the angular resolution MAGIC 
is as low as ~ 0.05°. This allows for a unprecedented 
resolution sky map from W51C. Figure [5] shows the 
MAGIC view of W51C overlapped with data from 
different wavelengths. The overall MAGIC emission 
above 700 GeV is contained within the region of 3 
counts above 1 GeV by Fermi/LAT (pink line in the 
figure). The higher angular resolution resolves that 
the gross of the VHE gamma-ray emission above 
700 GeV is spatially coincident with the shocked 
gas region reported by Koo et al. in [11[. This also 
supports the hadronic origin of the HE-VHE gamma 
rays. Furthermore there is an extension roughly 
following the shape of the remnants shell in radio and 
the Fermi/LAT HE emission. In addition this VHE 
emission towards the South-East goes in the direction 
of the PWN candidate. 
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Figure 5: Map of W51C region in different wavelengths. Shown in colors the MAGIC relative flux above 700 GeV 
(smooth by a Gaussian kernel of 0.065°) overlapped with the test statistics significance contours in black. The pink line 
shows the approximate contour of the HE emission detected by Fermi/LAT and the white dashed line shows the 
approximate contour of the H.E.S.S. VHE emission. The green dashed line shows the approximate contour of the SNR 
W51C. The dotted dark blue line shows the shocked gas region defined by Koo et al. [ll|]. On the left map the 
molecular clouds measured with the 66 km/s 13-CO line are shown in light blue. On the right map, the green contours 
show the radio data from Koo & Moon [l(J] in green. The X-ray source and PWN candidate CXO 192318.5+ 140335 and 
a 1720MHz OH maser are also shown in the map. 
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